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Abstract. We perform numerical hydrodynamic modeling of various physical processes that can form an HI ring 
as is observed in Holmberg I (Ho I). Three energetic mechanisms are considered: multiple supernova explosions 
(SNe), a hypernova explosion associated with a gamma ray burst (GRB), and the vertical impact of a high 
velocity cloud (HVC). The total released energy has an upper limit of ~ 10 54 ergs. We find that multiple SNe are 
in general more effective in producing shells that break out of the disk than a hypernova explosion of the same 
total energy. As a consequence, multiple SNe form rings with a high ring-to-center contrast K, <C 100 in the HI 
column density, whereas single hypernova explosions form rings with K, < 10. Only multiple SNe can reproduce 
both the size (diameter ~ 1.7 kpc) and the ring-to-center contrast (/C ~ 15 — 20) of the HI ring in Ho I. High 
velocity clouds create HI rings that are much smaller in size (<J 0.8 kpc) and contrast (/C <J 4.5) than seen in 
Ho I. We construct model position-velocity (pV) diagrams and find that they can be used to distinguish among 
different HI ring formation mechanisms. The observed pV-diagrams of Ho I (Ott et al. [20011 are best reproduced 
by multiple SNe. We conclude that the giant HI ring in Ho I is most probably formed by multiple SNe. We also 
find that the appearance of the SNe-driven shell in the integrated HI image depends on the inclination angle 
of the galaxy. In nearly face-on galaxies, the integrated HI image shows a ring of roughly constant HI column 
density surrounding a deep central depression, whereas in considerably inclined galaxies (i > 45°) the HI image 
is characterized by two kidney-shaped density enhancements and a mild central depression. 
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1. Introduction & Brinks l2UUT)l and Holmberg II (shell 21, according to 

Puche et al. lP39"2|l . 

There are a few nearby low-mass dwarf irregular galaxies 

(dlrr's) such as Holmberg I (Tully et al. 119781 Ott et al. There are a number of scenarios that could create 

EUrTT| . M81 dwA (Sargent et al. [P3g31 Westnfahl fc Puche HI holes in galactic disks (see Sanchez-Salcedo 2002), 

1993), Sagittarius DIG (Young & Lo I1997p . Sextans A which can loosely be divided into two groups compris- 

(Skillman et al. 119881 Stewart |l998), the HI morphology ing the energetic and non-energetic mechanisms. The en- 

of which is totally dominated by a single ring structure ergetic mechanisms assume a deposition of vast amount 

of size comparable to or bigger than their optical extent, of energy into the interstellar medium and include mul- 

The integrated HI image of these dlrr's shows a central tiple supernova explosions (Mac Low & McCray 119881 

HI hole surrounded by a denser HI ring. The contrast in De Young & Heckman 119941 Mac Low & Ferrara 119991 

the azimuthally averaged HI column density between the Silich et al.[2001), impact of high velocity clouds (Tenorio- 

central hole and the ring varies from ~ 3 - 4 in Sextans A Tagle et al. 119871 Comeron k Torra 11993 Rand & Stone 

and M81 dwA to ~ 15-20 in Holmberg I. HI rings that do I199rj[) . and gamma ray bursts (Efremov et al. [TS98 Loeb 

not dominate the overall HI extent of a galaxy, yet occupy & PernaEHEIl- The non-energetic mechanisms include a 

a considerable region (kpc scale), have also been found in combined action of thermal and gravitational instabili- 

other dlrr's such as DDO 47 (shell 13, according to Walter ties in the gas disk (Wada et al. 12000(1 . turbulent clearing 

(ElmegreenHinZI Walter & Brinks HM^ . and ultra-violet 
erosion of the HI disk (Vorobyov & Shchckinov I2()()4)l . It 



Send offprint requests to: E. I. Vorobyov appears however very unlikely that a central HI depres- 



E. I. Vorobyov, S. Basu: Formation of giant HI rings 



2 

sion surrounded by a denser HI ring with size compara- 
ble to the galaxy's optical extent could be produced by a 
non-energetic mechanism. Indeed, Vorobyov et al. (2004) 
have recently shown that the HI ring-like morphology of 
Holmberg I (Ho I) can be produced by multiple supernova 
explosions (SNe). 

In this paper, we consider two other energetic mech- 
anisms of HI ring formation, namely the impact of high 
velocity clouds (HVCs) and gamma ray bursts (GRBs), 
and show that they cannot explain the formation of the 
observed giant HI ring in Ho I. We generate the position- 
velocity diagrams and show that they can be used to dis- 
tinguish between the rings created by different energetic 
mechanisms. We find that the appearance of HI rings cre- 
ated by multiple SNe is sensitive to the inclination angle 
of the galaxy. In considerably inclined galaxies, the HI 
ring would rather appear as two kidney-shaped density 
enhancements similar to those observed in Sectans A. 

The paper is organized as follows. In Sect. El the nu- 
merical hydrodynamic model for simulating multiple SNe, 
GRBs, and the impact of HVCs is formulated. A compar- 
ative study of multiple SNe and GRBs of the same total 
energy is performed in Sect.0] The collision of HVCs with 
the galactic gas disk is considered in Sect. [31 The main re- 
sults are summarized in Sect.0 

2. Numerical model 

Our model galaxy consists of a rotating gas disk, stellar 
disk, and a spherically symmetric dark matter halo. The 
density profile of the stellar component is chosen as: 

p = p s0 sech 2 (z/,z s ) exp(-r/r s ), (1) 

where p s0 is the stellar density in the center of the galaxy, 
and z s and r s are the vertical scale height and radial scale 
length of the stellar component, respectively. We assume 
that the density profile of the dark matter (DM) halo can 
be approximated by a modified isothermal sphere (Binney 
& Tremaine[THE3) 



where the central density pho an d the characteristic scale 
length rh were given by Mac Low & Ferrara (J32JH2) and 
Silich & Tenorio-Tagle (2001) based on the study of 
the dark-to-visible mass ratios by Persic, Salucci, & Stel 

mm -. 

Pho = 6.3 x 10 10 V3 h-^ M kpc' 3 , (3) 

1 /2 

r h = 0.89 x 10~ 5 ( ^] h 1 ' 2 kpc. (4) 
V M ©/ 

Here, h is the Hubble constant in units of 
100 km s _1 Mpc -1 and Mh is the total halo mass. 
We adopt h — 0.65 throughout the paper. 



Our model galaxy is meant to represent Ho I and we 
use the observations of Ott et al. ( 2001) to constrain the 
parameters of our model. We note that Mh in Eqs. J3J 
and (0J) is, in fact, the total mass of the halo. We vary 
Mh until the actual halo mass confined within the HI di- 
ameter of Ho I (5.8 kpc, Ott et al. 2001J agrees with the 
observed value, 3.1 x 10 8 Mq. This fit yields a value of 
Mh = 6.0 x 10 9 Mq for Ho I. Once the total halo mass 
is fixed, we derive the parameters of the halo density dis- 
tribution using Eqs. |(3J) and Q, which are further used 
to compute the gravitational potential of the halo as de- 
scribed in Vorobyov et al. H2004[) . We have adopted a value 
of z s = 300 pc, which is typical for dwarf irregular galax- 
ies. The radial scale length r s — 1.7 kpc of the stellar 
disk is estimated from the 7 c -band radial surface bright- 
ness profile of Ho I. With r s and z s being fixed, the stellar 
density p° in Eq. ( is varied so as to obtain the mea- 
sured luminous stellar mass of M s = 1.0 x 10 s Mq. This 
results in p° w 0.02 Mq pc~ 3 . The adopted parameters 
of the stellar disk are further used to compute the stellar 
gravitational potential by solving the Poisson equation. 

Once the stellar and DM halo gravitational potentials 
are fixed, we obtain the initial gas density distribution by 
solving the steady-state momentum equation as described 
in Vorobyov et al. I|2004[l . We set the gas velocity disper- 
sion a = (RT/p) 1 / 2 to be 9 km s" 1 (Ott et al.EEU. We 
vary the rotation curve until the initial gas surface density 
distribution becomes exponential, which is in agreement 
with observations of the HI radial distribution in many 
dlrr's (Taylor et al. (1994). The observed rotation curve 
(RC) of Ho I is known to the accuracy of the inclination 
angle, which in turn is poorly determined (Ott et al.UOOl ) 
due to Ho Fs small inclination. We have chosen the initial 
RC in our model so that the initial radial gas surface den- 
sity distribution well reproduces the observed profile at 
r > 1 kpc, i.e. at radii which are not affected by the sub- 
sequent ring expansion. We do not expect our initial RC 
to match the currently observed RC, since the latter is al- 
ready affected by ring formation. Furthermore, accurately 
modeling the observed RC (for an assumed inclination an- 
gle) requires non-axisymmetric simulations, since the dy- 
namical center of Ho I is ~ 0.7 kpc offset from the morpho- 
logical center of the ring. This offset may complicate the 
shape of the observed RC, because it is measured around 
the dynamical center. This task is beyond the scope of 
our paper. The resulting radial gas surface density profile, 
as well as the Gaussian vertical scale height h of the gas 
distribution and initial gas rotation curve, are plotted in 
Fig- HI by the solid, dotted, and dashed lines, respectively. 
The total gas mass within the computational domain is 
M gas = 1.0 x 10 s M , of which ~ 30% is contributed 
by He (M m+Hc » 1.4 Mm, M = 1.27, Brinks ITM)|l . This 
value of M gas roughly agrees with that obtained by Ott et 
al. PM|| . 

A usual set of hydrodynamical equations in cylindrical 
coordinates (with the assumption of axial symmetry) is 
solved using the method of finite-differences with a time- 
explicit, operator-split solution procedure as used in the 
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Fig. 1. The equilibrium radial distribution of gas surface 
density (the solid line), the initial rotation curve (the 
dashed line) , and the Gaussian vertical scale height of gas 
distribution (the dotted line). 



ZEUS-2D code described in detail in Stone & Norman 
(1992J). The computational domain spans the range of 
7.2 kpc x 2 kpc in the vertical and horizontal directions, re- 
spectively, with a resolution of 5 pc. We have implemented 
the optically thin cooling curve given in Wada & Norman 
( 2001J) for a metallicity of one tenth of solar, which is typi- 
cal for dlrr's. The cooling processes taken into account are: 
(1) recombination of H, He, C, O, N, Si, and Fe; (2) col- 
lisional excitation of HI, CI-IV, and OI-IV; (3) hydrogen 
and helium bremsstrahlung; (4) vibrational and rotational 
excitation of H2; (5) atomic and molecular cooling due to 
fine-structure emission of C, C+, and O, and rotational 
line emission of CO and H2. We use an empirical heating 
function tuned to balance the cooling in the background 
atmosphere so that it maintains the gas in hydrostatic 
and thermal equilibrium; it may be thought of as a crude 
model for the stellar energy input. However, heating is pro- 
hibited at T gas > 2.0 x 10 4 K to avoid the effects of spuri- 
ous heating of a bubble's interior by the time-independent 
heating function. This is physically justified since most of 
the heating in the warm interstellar medium comes from 
the photoelectric heating of polycyclic aromatic hydrocar- 
bon molecules (PAHs) and small grains, which will be ei- 
ther evaporated or highly ionized in the bubble filled with 
hot supernova ejecta. Cooling and heating are treated nu- 
merically at the end of the time integration step using 
an implicit update to the energy equation. The implicit 
equation for energy density is solved by Newton-Raphson 
iteration, supplemented by a bisection algorithm for occa- 
sional zones where the Newton-Raphson method does not 
converge. In order to monitor accuracy, the total change in 
the internal energy density in one time step is kept below 
15%. If this condition is not met, the time step is reduced 
and a solution is again sought. In the following we give 
a brief description on the mechanisms of energy injection 
used in our numerical simulation. 



3. Energetic mechanisms of HI ring formation 

Here, we briefly review possible mechanisms of HI ring 
formation as observed in Ho I. 

Consecutive supernova explosions. The origin of the gi- 
ant HI rings in dlrr's is traditionally thought to lie in the 
combined effect of stellar winds and supernova explosions 
produced by young stellar associations (see e.g. Ott et 
al. [2001). In our simulations, the energy of supernova ex- 
plosions is released in the form of thermal energy in the 
central region with a radius of 30 pc filled with the hot 
(T ~ 10 7 K) and rarified (n ~ 10~ 3 cm -3 ) gas, which is 
presumably formed by the previous action of stellar winds. 
We use a constant wind approximation, i.e. at each time 
step we add energy to the source region at a rate of E = L, 
where L is the mechanical luminosity defined as the total 
released energy of SNe divided by the duration of the en- 
ergy input phase. We choose the energy input phase to 
last for 30 Myr, which corresponds roughly to the differ- 
ence in the lifetimes of the most and least massive stars 
capable of producing SNe in a cluster of simultaneously 
born stars. Since in the present simulations we deal with 
large stellar clusters with hundreds of supernovae, the re- 
lease of the energy of SNe in the form of thermal energy 
is justified (Mac Low & McCrav lT^jl . 

Gamma ray bursts. The consecutive SN explosions may 
not be the only mechanism that could release 10 53 — 
10 54 ergs of energy, enough to form giant HI rings. Another 
mechanism has been suggested by Efremov et al. (1998) 
and Loeb & Perna l|1998[) . wno argued that the GRB 
explosions are powerful enough to make kpc-size shells 
in the interstellar media of spiral and irregular galaxies. 
Although the physics of GRBs is still poorly understood, 
the general picture emerging is that they are highly en- 
ergetic events (< 10 54 ergs, see e.g. Paczvhski 119981 that 
release energy in a short period of time (of the order of 
a few seconds). Consequently, we model the GRB explo- 
sion (hereafter, "hypernova" according to Paczyhski 1998) 
by an instantaneous release of thermal energy within a 
sphere of 30 pc, filled with hot (T ~ 10 7 K) and rarified 
(n ~ 10~ 3 cm~ 3 ) gas. We have also explored the injection 
of hypernova energy in the form of kinetic energy and 
found that it does not noticeably influence the dynamics 
of the hypernova-driven shell. 

Impact of high velocity clouds. Another possible mech- 
anism that could create HI holes in the galactic disks was 
proposed by Tenorio-Taglc et al. (1987), who argued that 
the infall of high velocity clouds could deposit 10 52 to 10 54 
ergs per collision. The numerical simulations of Tcnorio- 
Tagle et al. Ijl987|) have indeed shown that HVCs are ca- 
pable of forming the giant curved arcs and cavities in the 
Galactic disk (see also Rand & Stone (1990.) for numerical 
simulations of HVC impact in NGC 4631). Most previous 
numerical simulations have been concerned with the in- 
fall of HVCs in massive galaxies, because their collisional 
cross section is much larger than that of dwarf irregulars. 
Moreover, the gas disks in dwarf irregulars are in general 
thicker than those of massive spiral galaxies, which would 



4 



E. I. Vorobyov, S. Basu: Formation of giant HI rings 





v 


\ 








1=18 Myr 




t=27 Myr 


t=4S Myr 


t=72 Myr 


t=1 10 Myr 



Galactocentric radius (kpc) 

Fig. 2. Temporal evolution of the gas volume density dis- 
tribution after the release of 2 x 10 53 ergs of thermal energy 
at z — 0, r = pc. The upper panels correspond to 200 
consecutive SN explosions over 30 Myr, while the lower 
panels show the impact of a single hypernova of the same 
total energy of 2.0 x 10 53 ergs. The scale bar is in gm cm~ 3 . 

make it more difficult for an HVC to penetrate the disk. 
In summary, the formation of a giant HI hole surrounded 
by a denser HI ring (as is observed in Ho I and other 
dwarf irregulars) by the infall of HVCs is not obvious and 
requires further investigation. 

Taking into account the diameter of the HI ring in 
Ho I (~ 1.7 kpc), we consider the most energetic HVCs, 
which cover a velocity range 200 km s~ x ■ 

W HVC 1_ 

300 km s _1 and have HI column density 10 20 cm~ 3 < 
Nm < 10 21 cm~ 3 . In our simulations the kinetic energies 
of HVCs range from 1.0 x 10 53 to 2.5 x 10 54 ergs, which 
corresponds to a variation in HVC masses of 2.5 x 10 5 Mq 
to 2.5 x 10 6 Mq. Such massive and energetic HVCs most 
probably have an extragalactic origin. 
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Fig. 3. The same as Fig. [21 but for the energy input of 
5 x 10 53 ergs. 
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Fig. 4. The same as Fig. [21 but for the energy input of 
10 54 ergs. 



4. Supernova explosions versus gamma ray bursts 

We start by showing in Fig. [5] the temporal evolution of 
the distribution of the gas volume density produced by 200 
consecutive SNe (the upper panels) and by a hypernova 
explosion of the same total energy of 2 x 10 53 ergs (the 
lower panels). The overall gas dynamics is similar in both 
cases. However, there are minor differences seen not only 
in the dynamics of the hot gas but also in the shape of 
the blown-up shell. In the case of a hypernova explosion, 
there is no "blowout" observed (i.e. the shell breaking out 
of the disk and pumping the hot gas into the intergalactic 
medium) . The hot gas is always confined inside the shell in 
the early expansion phase and it cools down at later times 
t > 40 Myr due to the radiative cooling. As a result, the 
shell collapses at t — 45 Myr, creating a mildly compressed 
central core, the gas density of which is however below the 
Jeans limit. In the case of multiple SNe, part of the hot gas 
is lifted to a higher altitude of z > 1 kpc due to a buoyancy 
effect developing at t ~ 45 Myr. It is also seen that in the 



early expansion phase the shell is on average thicker in 
the case of the hypernova explosion than in the case of 
multiple SN explosions. Since the numerical modeling of 
Vorobyov et al. (J2004 ) suggested that the shell in Ho I has 
already blown out of the disk, we do not consider further 
the energy release of 2.0 x 10 53 ergs. 

The difference in the shell dynamics between the multi- 
ple SN and hypernova becomes pronounced as one consid- 
ers more energetic explosions. For instance, Fig. shows 
the distribution of the gas volume density produced by 
500 consecutive SN explosions (the upper panels) and by 
a single hypernova explosion of the same total energy 
5.0 x 10 53 ergs (the lower panels). Now, the dynamics of the 
SN-driven shell shows a clear blowout phase at t « 45 Myr 
after the beginning of the energy input phase. On the con- 
trary, the hypernova-driven shell never breaks out of the 
disk, though its energy is equal to the total energy released 
by multiple supernovae. This difference remains for pre- 
sumably an upper limit energy release of 10 54 ergs that 
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Fig. 5. The vertical cut through the distribution of the 
gas volume density produced by 500 consecutive SNe at 
t=27 Myr (the left panel) and single hypernova of the 
same total energy of 5.0 x 10 53 ergs (the right panel). The 
gas velocity field is shown by the arrows. 



a stellar cluster could produce in Holmberg I. As seen 
in the upper panels of Fig. ^ the shell dynamics gov- 
erned by 1000 consecutive SN explosions exhibits a vio- 
lent blowout, whereas that governed by a hypernova of 
the same total energy of 10 54 ergs shows almost no sign 
of blowout. This is in good agreement with the previous 
model of Efremov, Ehlerova & Palous (1999), who used a 
thin shell approximation to study the evolution of shells 
formed by a single hypernova and multiple SNe. These au- 
thors have also found that an abrupt energy input creates 
shells that do not blow out to the galactic halo for energies 
<J 10 54 ergs. Noticeably, the final fate of the gas distribu- 
tion is similar in all cases: the hot bubble fills in and the 
gas disk mostly recovers its pre-explosion appearance af- 
ter > 160 Myr with a slightly more centrally condensed 
radial gas distribution. 

The occurrence of the blowout in the case of multiple 
supernova explosions and its absence in the case of a hy- 
pernova explosion can be understood if one considers the 
dynamics of the hot gas filling the shell interior. In Fig. 
we plot the gas volume density distribution and veloc- 
ity field produced by 500 consecutive SN explosions (the 
left panel) and a hypernova of the same total energy of 
5.0 x 10 53 ergs (the right panel). An obvious difference is 
seen: the hot gas ejected by multiple SNe forms a "vortex" 
that acts to de-stabilize the swept-up shell of cold mate- 



rial via a Kelvin-Hclmholtz instability. The shell loses its 
smooth elliptical form and develops "ripples" at its tops. 
The dynamical pressure of the hot gas (note that its ve- 
locity is much higher than in the case of a hypernova ex- 
plosion) accelerates the "rippled" shell. As a consequence, 
the shell shows a strong Rayleigh- Taylor instability via the 
development of a characteristic 'spike-and-bubble' mor- 
phology, as is seen in the left panel of Fig. [5] In the case 
of a hypernova explosion of the same total energy, the 
gas velocity field plotted in the right panel of Fig. lacks 
any circular motion inside the shell. The smooth elliptical 
shell of cold material of roughly the same size as in the 
case of multiple SNe expands at roughly the sound speed. 
The shell has already started losing its pressure support 
(as implied by the thickening of the shell walls) due to 
radiative cooling of its hot interior. As a consequence, the 
decelerating shell does not develop a Rayleigh- Taylor in- 
stability and the blowout phase is not observed. We note 
here that a hypernova-driven shell expands out to a certain 
radius faster than a shell created by multiple SNe with the 
same total energy, as reported by Efremov et al. (^999). 
Our numerical simulations have shown that the difference 
in the dynamics of shells produced by multiple SNe and 
a single hypernova of the same total energy intensifies as 
one considers more flattened gas systems. The maximum 
difference is found in the strongly vertically stratified gas 
distributions assumed for massive spiral galaxies. 

The vortex forms since the shell becomes elongated in 
the vertical direction. Consider the simplified situation of 
hot gas ejected radially by SNe which is reflected from 
the cold dense walls of the elliptical shell. The gas will 
not be reflected along the local normal direction to the 
shell unless it is moving along z = or r = 0. In gen- 
eral, a tangential component of the velocity of hot gas Vt 
(with respect to the shell's walls) is generated. It is always 
directed upwards for gas above the midplane and down- 
wards for gas below the midplane. Let us consider the 
upper hemisphere. The occurrence of vt ^ near the shell 
walls and the axial symmetry of the shell makes the hot 
gas (streaming upwards along the shell walls) accumulate 
at the top of the shell, because it cannot pass through the 
symmetry axis. Both the growing pressure of this hot gas 
(due to its negligible cooling) and the downward pull of 
the z-componcnt of the galactic gravitational field lead to 
a downward flow along the axis of symmetry. This com- 
pletes the circle and generates a vortex structure as seen in 
the left panel of Fig. The development of such vortices is 
also seen in the axially symmetric numerical simulations of 
Recchi et al. (2001). However, a strong non-axisymmetry 
of the SNe-driven shell may complicate the formation of 
the vortex. 

The driving force for the vortex is a continuous or 
quasi-continuous release of energy by SNe of a stellar clus- 
ter or a group of closely located stellar clusters. This is also 
the reason why vortices do not form in a shell created by a 
single hypernova explosion. We have found that the vortex 
develops even when the number of SNe is quite moderate 
(~ 20) and the release of energy is discrete. The length 
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scale of the vortex is approximately equal to the shell's 
semi-minor axis and time scale is limited by the duration 
of energy input from SNe, i.e. < 30 Myr. In general, at 
least two vortices can co-exist within a single shell: one 
in the upper hemisphere and the other in the lower hemi- 
sphere of the shell. 

Hypernovae tend to form complete shells. Such shells 
would appear in nearly face-on dwarf galaxies as HI rings 
with a low ring-to-center contrast in the HI column den- 
sity. For instance, in Fig. we plot the contrast in the 
HI column density between the ring and the central de- 
pression, JC = N r (HI)/N c (HI) 1 as a function of the ring 
radius R obtained for three different energy inputs. The 
ring column density N r (HI) is computed by azimuthally 
averaging N(HI) around the ring. An inclination angle 
(i.e. the angle between the rotation axis of the galaxy and 
the line of sight) of i = 5° is assumed, which is appropriate 
for a nearly face-on galaxy. The solid lines give the con- 
trast JC for the multiple SN explosions, whereas the dashed 
lines do that for the hypernova explosion. As is seen, the 
ring-to-center contrast JC produced by the hypernova ex- 
plosion never exceeds 10 even for the upper limit energy 
of hypernova explosions of 10 54 ergs. On the other hand, 
multiple SNe can form HI rings with a much higher ring- 
to-center contrast JC <; 100. The latter takes place when 
the SN-driven shell is in the blowout phase. Note that the 
difference in JC between the HI rings produced by multiple 
SNe and those produced by hypernovae smears out if the 
total released energy is < 2 x 10 53 ergs. 

4.1. PV-diagrams 

The position-velocity (pV) diagram is a powerful tool to 
search for coherent structures in the interstellar medium. 
In this section we compute model pV-diagrams in order 
to determine their utility for distinguishing between the 
SN-driven shells and those produced by a hypernova ex- 
plosion of the same total energy. The HI flux density Shi is 
obtained using the following conversion formula (Binney 
& Merrifieldf 



Mi 



HI 



A; 



2.35 x 10 5 



( D 



Shi 



Mr 



km f 



(5) 



VMpc 

where for the distance D we take 3.5 Mpc, appropriate for 
Holmberg I. The position- velocity cuts are taken along the 
ma jor axis of the projected shell. The width of the cuts is 
30 pc. We tried twice wider cuts (60 pc), but found little 
difference in the appearance of the resulting pV-diagrams 
as compared to those obtained with narrower cuts. The 
gas is assumed to be thermalized when constructing the 
model pV-diagrams, i.e. a Maxwellian velocity distribu- 
tion is assumed for the gas in each computational cell, 
with the scale determined by the local gas temperature 
(see Mashchenko & Silich 1995). No attempt was made 
to include the effects of turbulence. The model galaxy is 
assumed to have an inclination angle of 5°. The position- 
velocity cuts along the minor axis of the projected shell 
have similar appearance for such low inclination angles. 
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0.6 0.8 
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Fig. 6. The contrast in the HI column density between 
the ring and the central depression JC = N T (HI)/N C (HI) 
as a function of the ring radius. The solid lines show the 
contrast JC produced by multiple SNe, whereas the dashed 
lines show that produced by a single hypernova of the same 
total energy. An inclination angle of 5° is assumed. The 
energy input is indicated in the right upper corner of each 
panel. The dashed-dotted and dotted lines in the lower 
panel give the contrast JC for an inclination angle of 45°. 



Figurc[7|shows the model pV-diagrams of the shell pro- 
duced by 500 consecutive SNe. The quantity Shi is plotted 
in all model pV-diagrams. Two bright blobs elongated in 
the vertical direction are apparent in each pV-diagram. 
They represent the dense walls of the shell expanding in 
the plane of the galaxy, i.e. perpendicular to the line of 
sight at the adopted inclination angle of 5°. If the shell 
has not yet broken out of the disk (t = 9 and 18 Myr), 
these two blobs appear to close at their tops and form 
a complete elliptical ring. The relative amplitude of Shi 
around the ring is however varying by roughly an order 
of magnitude, reflecting the difference in the HI column 
density along the major axis of the projected shell - there 
is much less gas on the approaching and receding sides of 
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Fig. 7. The pV-diagrams of the shell created by 500 con- 
secutive SNe. The pV-cuts are taken along the major axis 
of the model galaxy viewed at an angle of 5° . The evolu- 
tionary time of the shell is indicated in the right upper cor- 
ner. The quantity .Shi is plotted in all model pV-diagrams. 

the shell than on those expanding perpendicular to the 
line of sight . 

The elliptical ring is virtually absent in the later ex- 
pansion phase at t — 36 and 45 Myr when the shell has 
already broken out of the disk (the lower frames in Fig.0). 
Instead, the pV-diagrams are totally dominated by two 
high-intensity blobs. A few smaller patches seen between 
the blobs represent the dense, cold clouds fragmented off 
the shell and pushed by the dynamical pressure of the hot 
gas to a higher altitude of > 1 kpc. The thickening of 
the shell after the blowout is also noticeable in the lower 
frames of Fig. The same tendency is found in the pV- 
diagrams of the shell created by 1000 consecutive SNe. 
The complete ring is seen only in the very early phase of 
the shell expansion at t=9 Myr and is virtually absent 
after the blowout at t = 27, 45, and 54 Myr. 

Figure [S] shows the pV-diagrams of the shell produced 
by a hypernova of total energy 5 x 10 53 ergs. The position- 
velocity cuts are taken along the major axis of the pro- 
jected shell at four different evolutionary times: t = 9, 
18, 27, and 36 Myr, i.e., before it has expanded out to 
its maximum size (t = 9 and 18 Myr) and after it has 
started to collapse due to the radiative cooling of the shell 
interior. An inclination angle i = 5° is assumed. The ellip- 
tical hole surrounded by a ring of higher HI flux density 
Shi is clearly visible in those pV-diagrams. The HI flux 
Shi from the ring is ~ 7 — 10 \i Jy, which should, in prin- 
ciple, be detectable. Hypernova explosions tend to form 
complete shells that are clearly visible in the pV-diagrams. 
This characteristic feature of a hypernova-driven shell is 



also present when one considers very energetic explosions. 
For instance, Fig. shows the pV-diagrams of the shell 
produced by a hypernova of the total energy of 10 54 ergs. 
The elliptical hole surrounded by the ring is clearly seen 
even when the shell has already started to collapse at 
t = 45 Myr, as implied by the presence of a vertical bar 
at zero radial offset (see the lower right panel in Fig. |5J). 
Note the size of the shell in Figs. [3] and It is twice 
the local Gaussian scale height h, which is ~ 350 pc at a 
galactocentric radius of 0.5 kpc in our model galaxy. 

The comparison of Fig.[7|and Figs. 03 El indicates that 
the SNe-driven shell in its late expansion phase (when 
its radius exceeds 1.5 h) may be distinguished from that 
produced by the single hypernova explosion based on 
its appearance in the pV-diagrams. The pV-diagrams of 
the hypernova-driven shell show the characteristic ellipti- 
cal ring in virtually any evolutionary phase of the shell, 
whereas the pV-diagrams of the SN-driven shell do that 
only in the early expansion phase, well before the blowout. 
The present numerical simulations (see also Vorobyov & 
Shchekinov 2004) show that the SN-driven shell breaks 
out of the disk when its radius exceeds ~ 1.5 h . As a con- 
sequence of the blowout, the SN-driven shell transforms 
into an open cylinder and the pV-diagrams show two par- 
allel blobs representing the walls of the expanding cylin- 
drical shell. On the contrary, the hypernova-driven shell 
appears to survive the expansion even when its radius ex- 
ceeds 2.0 h. As a result, the characteristic elliptical ring 
representing the complete expanding shell is clearly seen 
in the pV-diagrams. Thus, if the pV-diagrams show an 
elliptical hole surrounded by a denser ring, the radius of 
which is greater than 1.5 h (as in the lower frames of Fig|Hl 
and El, this may indirectly point to a hypernova origin of 
the shell. 



4.2. Effects of inclination 

The model pV-diagrams considered in Sect. 14. ll were con- 
structed under the assumption that the model galaxy is 
viewed nearly face-on, namely at an inclination angle of 
5° . In this section we study how a higher inclination angle 
could influence the appearance of the SN- and hypernova- 
driven shells in the pV-diagrams. Figure EH shows the 
model pV-diagrams of the shell created by 1000 consecu- 
tive SNe. The pV-diagrams are constructed at three dif- 
ferent evolutionary times t — 9, 27, and 45 Myr after the 
beginning of the energy injection phase. The left and right 
panels in Fig. ^| are the position- velocity cuts along the 
minor and major axes of the projected shell, respectively. 
An inclination angle of i = 45° is adopted. The character- 
istic hole surrounded by a ring of higher HI flux density 
(see Sect. 14 .l|) is barely seen even in the very early expan- 
sion phase at t — 9 Myr due to the projection smearing at 
higher inclination angles. In the later expansion phase at 
t = 45 Myr the pV-cuts taken along the major axis show 
two elongated blobs, the morphological centers of which 
are symmetrically displaced by |Au| ^7 — 8 km s^ 1 with 
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Fig. 8. The pV-diagrams of the shell produced by a single 
hypernova of the total energy of 5.0 x 10 53 ergs. The pV- 
cuts are taken along the major axis of the model galaxy 
viewed at an inclination angle of 5°. 
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Fig. 9. The same as Fig. 00 but for a single hypernova of 
the total energy of 10 54 ergs. 



respect to v = km s _1 . The displacement is defined as 
Aw = ±u ro t sin(z), where v ro t is the rotation velocity of 
the SN-driven shell. The positive and negative Av corre- 
spond to the receding and approaching parts of the shell, 
respectively. It is important to note here that similar pV- 
diagrams are obtained for Ho I by Ott et al. 1)20011 see 



Fig. 10. The pV-diagrams of the shell created by 1000 
consecutive SNe. The pV-cuts are taken along the minor 
(the left panels) and the major (the right panels) axis of 
the model galaxy viewed at an angle of 45°. 



their Fig. 8). A large displacement between the blobs with 

is 



>, respect to the systemic velocity of v sys = 141.5 km s~ 
, a clearly seen in the pV-cut taken along the major axis of 
> Ho I. However, we found it difficult to estimate the value of 
Av from the observed pV-diagram shown in Fig. 8 of Ott 
et al. H20()l|l due to its considerable asymmetry (caused 
probably by an off-center location of the stellar cluster re- 
r sponsible for the creation of the shell). The modeled value 
of Av depends also on the adopted rotation curve of Ho I, 
which is poorly known (Ott et al. I2001fl . Therefore, it is 
problematic to obtain the inclination angle of Ho I by 
comparing the modeled and observed Av. Nevertheless, 
we conclude that Ho I is viewed at an inclination angle 
considerably higher than i = 5°, because the modeled pV- 
cuts of Fig. [7| obtained for an inclination angle of only 5° 
do not show a noticeable displacement. 

The pV-diagrams of the shell created by a single hy- 
pernova explosion of the total energy 10 54 ergs are shown 
in Fig. ^2 The left panels give the pV-cuts along the mi- 
nor axis of the projected shell, whereas the right panels 
do those along the major axis. The pV-diagrams of the 
shell created by a single hypernova explosion are notice- 
ably different from those created by multiple SNe. The 
hypernova-driven shell does not break out of the disk even 
when its radius exceeds 2.0 h. As a consequence, the shell 
preserves its elliptical form and the characteristic elliptical 
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Fig. 11. The same as in Fig llOl but for a single hypernova 
of the same total energy of 10 54 ergs. The pV-cuts taken 
along the minor axis of the hypernova-driven shell show 
the characteristic S-shaped structure due to an expansion 
in the vertical direction. 



hole surrounded by the ring of a higher HI flux intensity 
is apparent in the pV-cuts taken along the major axis 
(the right panels in Fig. 111(1 at t = 9 and 27 Myr after 
the explosion. The elliptical ring dissolves at t > 45 Myr 
due to the projection effect (note that in nearly face-on 
galaxies the elliptical ring can still be resolved at that 
phase). In that later phase, the hypernova-driven shell 
can barely be distinguished from that created by multiple 
SNe based on the appearance in the pV-diagrams taken 
along the major axis of the projected shell. Fortunately, 
the pV-cuts taken along the minor axis of the hypernova- 
driven shell show the characteristic S-shaped structure. 
The pV-cuts along the minor axis are sensitive only to 
the expansion/contraction motion of the shell. As is seen 
in the lower panels of Fig. the hypernova-driven shell 
at t = 27 — 45 Myr expands in the vertical direction, but 
shows the lack of expansion or contraction motion in the 
galactic plane (except for the very late phase when the 
shell has already started to collapse). As a result, the shell 
seen in the pV-diagrams as the elliptical ring transforms 
into the S-shaped structure in the later expansion phase. 
Note that neither elliptical rings nor S-shaped structures 
are seen in the observed pV-diagrams of Ho I plotted in 
Fig. 8 of Ott et al. p¥)T)l . 



As was demonstrated in Fig. El the HI rings created 
by multiple SNe have a much higher contrast in the HI 
column density between the ring and central depression, 
K. = N r (HI)/N c (HI), than those created by a single hy- 
pernova explosion of the same total energy. However, a 
higher inclination may considerably reduce the value of 
K, due to the projection smearing effect. For instance, the 
dotted-dashed line in Fig. HJ shows the contrast JC as a 
function of the ring radius R produced by 1000 consec- 
utive SNe, whereas the dotted line does that for the hy- 
pernova of the same total energy. An inclination angle of 
45° is assumed. The ring-to-center contrast in N(HI) pro- 
duced by both multiple SNe and hypernova stays around 
JC = 2 — 3. The only exception is seen when the SN-driven 
shell expands out to R ss 0.7 kpc and breaks out of the 
disk. Then, the ring-to-center contrast reaches a maximum 
value of JC w 8. It is instructive to compare the ring-to- 
center contrast and size of HI rings produced in our simu- 
lations (by both multiple SNe and single hypernovae) with 
the observed contrast JC and size of the HI ring in Ho I. 
Fig. 1121 shows the observed HI map of Ho I as obtained by 
Ott et al. (200! i. The diameter of the ring is ~ 1.7 kpc 
and the ring-to-center contrast is JC ~ 15 — 20. Our simula- 
tions show (see Fig. |SJ) that both multiple SNe and single 
hypernovae may form rings of up to 2 kpc in diameter. 
However, single hypernovae fail to produce rings with the 
observed contrast of JC ~ 15 — 20. This is true even in the 
most favorable case of near zero inclination angle. Taking 
into account that both the observations of Ott et al. and 
numerical modeling of Vorobyov et al. (|2004[1 set the lower 
limit on the inclination angle of Ho I as % < 15°, the forma- 
tion of the HI ring with such a high contrast JC ~ 15 — 20 
by means of a hypernova explosion seems highly problem- 
atic. Apparently, only multiple SNe can produce HI rings 
with such a high ring-to-center contrast JC. 

The integrated HI image of the shell created by mul- 
tiple SNe appears in nearly face-on galaxies as a ring of 
elevated HI density surrounding the central depression. 
For instance, the left upper panel in Fig. ^| shows the 
integrated HI image of the shell created by 500 consec- 
utive SNe. The HI image is taken at t = 45 Myr after 
the beginning of the energy input phase, when the shell 
has already broken out of the disk. An inclination angle 
of i — 5° is assumed. As is obviously seen, the blown-out 
shell in nearly face-on galaxies appears as a ring with a 
high ring-to-center HI column density contrast (~ 25—30). 
There is little azimuthal variation in the HI column den- 
sity around the ring. However, the integrated HI image 
of the blownout shell observed in substantially inclined 
galaxies shows a considerable variation in the HI column 
density around the ring. For instance, at higher inclination 
angles of 45° and 60° the relative amplitude of azimuthal 
variations in A (HI) around the ring is k 2 and w 3, re- 
spectively. In this case, the integrated HI image is dom- 
inated by two kidney-shaped density enhancements and 
shows a mild central minimum. It is interesting to note 
that the Sextans A and M81 dwA dwarf irregular galaxies 
have a somewhat similar HI distribution (see Skillman et 
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Fig. 12. Integrated HI emission of Ho I. In the bottom 
left corner we show the half-power beam (8'.'2 x 7. (/.'). The 
grayscale represents the HI column density. The data have 
been obtained from the Very Large Array observations 
(Ott et al.EHHIl. 

al. 119881 Stewart [1998 ). In the case of Sextans A, the rel- 
ative amplitude of azimuthal variations in the HI column 
density around the ring is « 2.5 — 3 (according to Skillman 
et al. I1988J) and the inclination angle determined from the 
outer HI contours is i = 36° ± 4° (assuming the circular 
symmetry). Note that the narrow ring-shaped structures 
seen in Fig. 1131 are in fact the dense fragments of a broken- 
up shell, which appear as rings rather than dense clumps 
due to the axisymmetric nature of our numerical simula- 
tions. In edge-on galaxies, the integrated HI image would 
resemble a 'dumbbell', as shown in the right lower panel of 
Fig. 1131 The absence of noticeable azimuthal variations in 
HI column density around the HI ring in Fig. El suggests 
that Ho I is viewed at an inclination angle considerably 
smaller than i = 45°. This is in agreement with previous 
estimates by Ott et al. l|2001[l and Vorobyov et al. 1)2UU4[1 . 
who found that Ho I has an inclination angle of < 20°. 

4.3. Off-plane explosions 

A general assumption made in the above simulations of 
both SN-driven and hypernova-driven shells was that the 
energy is released exactly in the midplane of the model gas 
disk. However, taking into account a considerable thick- 
ness of the gas disks in dlrr's, the occurrence of a hy- 
pernova or massive stellar cluster at a position moder- 
ately offset with respect to the galactic midplane, appears 
plausible. In this section we investigate how the off-plane 
explosions may influence the appearance of the shells in 
the pV-diagrams. 

If the stellar cluster is located at a moderate height of 
z — 100 pc (~ 1/3 h) above the midplane of the galaxy, 
the temporal evolution of the gas volume density is found 
to be similar to that obtained in Sect. 01 for the midplane 
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Fig. 13. The integrated HI image of the ring produced 
by 500 consecutive SNe. The inclination angle is indicated 
in the right upper corner. The scale bar is in cm -2 . The 
origin of narrow rings is an assumption of axisymmetry in 
our numerical simulations. 

stellar cluster. For instance, the upper panels in Fig. ITU 
show five temporal snapshots of the gas volume density 
distribution created by the energy release of 1000 con- 
secutive SNe. A more intensive blowout takes place in the 
upper part of the galaxy where the stellar cluster is nested. 
However, the pV-diagrams and ring-to-center contrast in 
HI column density /C are only slightly affected by this mild 
asymmetry. 

The temporal evolution of the shell created by a sin- 
gle hypernova that is moderately offset with respect to 
the midplane of the galaxy produces some noticeable dif- 
ferences. The lower panels in Fig. 1141 give five temporal 
snapshots of the gas volume density distribution created 
by the hypernova of total energy 10 54 ergs that is placed 
100 pc above the midplane of the galactic gas disk. The 
egg-like shell is formed after « 15 Myr, the upper part of 
which retains a relatively smooth appearance for several 
tens of Myr before it breaks due to a Ray leigh- Taylor in- 
stability and dissolves at « 70—80 Myr after the explosion. 
Note that in the case of multiple SNe such smooth kpc- 
sized arcs of dense material are never formed in our nu- 
merical simulations. Instead, inspection of Figs. 121 II and 
Fig-Elindicates that multiple SNe tend to produce patchy 
filaments and dense clumps of a few hundred parsecs in 
size (except for the very late phase when the shell has al- 
ready collapsed). Indeed, a Rayleigh- Taylor instability, the 
effect of which is strengthened by the preceding Kelvin- 
Hclmholtz instability, efficiently fragments the SN-driven 
shell when the latter breaks out of the gas disk. Longer 
fragmentation times of hypernova-driven shells as com- 
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Fig. 14. Off-plane explosions. Temporal evolution of the 
gas volume density distribution after the release of 1.0 x 
10 54 ergs of thermal energy at r = 0, z = 100 pc. The up- 
per panels correspond to 1000 consecutive SN explosions 
over 30 Myr, while the lower panels show the impact of 
a single hypernova of the same total energy of 10 54 ergs. 
The scale bar is in gm cm" 3 . 



pared to fragmentation times of shells created by multiple 
SNe were also found in the model of Efremov et al. ([1999). 
It is interesting to mention in this context the W4 HII re- 
gion in the Perseus arm of the outer Galaxy. Dennison et 
al. I|1997[) have observed a highly elongated shell of Ha 
emission, the base of which is located near the W4 HII re- 
gion. The Ha shell appears to close far above the Galactic 
plane and is interpreted by Basu et al. Q1999[) to be the 



dense ionized wall of swept-up gas surrounding a bubble 
of hot rarefied gas created by the stellar winds from the 
massive stars. However, the apparent smooth appearance 
of the Ha shell implies that the Rayleigh- Taylor insta- 
bility is suppressed there, possibly due to the stabilizing 
effect of a swept-up tangential magnetic field in the shell 
(Komljenovic, Basu, & Johnstone 1999). An alternative 
explanation may involve a shell formed by a hypernova 
explosion. In that case, the characteristic growth time of a 
Rayleigh- Taylor instability in the expanding shell is much 
longer than in the case of multiple SNe, and the shell may 
survive for a longer time and expand to a higher altitude. 

Off-plane hypernova explosions may leave noticeable 
signatures in the pV-diagrams. While the upper part of 
the egg-like shell is thin and contains a small fraction of 
the total gas mass, the lower part of the shell accumulates 
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Fig. 15. Off-plane explosion. The pV-diagrams of the 
shell created by a single hypernova of the total energy 
of 10 54 ergs. The pV-cuts are taken along the major axis 
of the model galaxy viewed at an inclination angle of 5°. 



a larger amount of gas mass and becomes clearly visible in 
the pV-diagrams. FigureEDshows the pV-cuts taken along 
the major axis of the shell created by a single hypernova 
of the total energy of 10 54 ergs, which is placed 100 pc 
above the midplane of the galaxy. An inclination angle of 
5° is adopted. The bottom of the egg-like shell appears in 
the pV-diagrams as a bright arc connecting two vertical 
blobs. The latter are the walls of the shell expanding in the 
horizontal direction. The top of the shell is clearly visible 
only in the very early expansion phase at t < 10 Myr. 
However, it could be seen in Ha emission, if a powerful 
ionizing source is located inside the shell, as it might be 
the case for W4 HII region discussed above. 



5. High velocity clouds 

In this section we numerically study the vertical impact 
of high velocity clouds as a potential mechanism of the 
giant HI ring formation in Ho I. The numerical simula- 
tions of multiple SNe in Sect. 01 have shown that it re- 
quires > 2.0 x 10 53 ergs of released energy to form HI 
rings with HI ring-to-center contrast JC > 10. Hence, 
we consider only very energetic HVCs that could poten- 
tially deposit a comparable amount of kinetic energy to 
the interstellar medium. We have performed a parame- 
ter study assuming that HVCs cover a range of veloc- 
ities 200 km s" 1 < vhvc < 300 km s" 1 and masses 
2.5 x 10 5 M Q < M H vc < 2.5 x 10 6 M Q . Thus, the kinetic 
energies of HVCs range from 1.0 x 10 53 to 2.5 x 10 54 ergs. 
Clouds are assumed to have a constant density and be ei- 
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Fig. 16. The distribution of the gas volume density cre- 
ated by the collision of a cloud of Mhvc = 2.5 x 10 6 M© Fig. 17. The integrated HI image of the gas volume den- 
and whvc = 200 km s _1 with the galactic gas disk. sity distribution shown in Fig.EO An inclination angle of 

5° is adopted. The scale bar is in cm~ 2 . 



ther spherical or cylindrical in shape. The initial position 
of the cloud center is 3.2 kpc above the midplane. 

Figurc^j]shows the temporal evolution of the distribu- 
tion of the gas volume density produced by the impact of 
a cloud of Mhvc = 2.5 x 10 6 M Q and w H vc = 200 km s _1 . 
The total kinetic energy of the cloud is thus equal to 
1.0 x 10 54 ergs. The maximum column density iV(HI) of 
the cloud measured along its diameter is 4.3 x 10 20 cm~ 2 . 
The impinging HVC creates a hot rarified bubble by push- 
ing matter toward the galactic plane. The bubble does not 
form a complete sphere and fills out ~ 40 Myr after the 
onset of the collision. The HVC continues to push the gas 
along its way and eventually breaks through the disk. At 
this stage (t w 110 Myr), the initial shape of the HVC 
is completely lost due to a combined action of Rayleigh- 
Taylor and thermal instabilities. The HVC emerges above 
the galactic plane as a complicated network of dense fila- 
ments, clumps, and arcs. The formation of dense filaments 
and clumps was also reported in numerical simulations by 
Franco et al. ( 1985J). It takes another 70 Myr for this rem- 
nant to fall back onto the disk and create another bubble, 
though less prominent than the first one. 

Although the kinetic energy of the impinging cloud is 
high (10 54 ergs) and the cloud produces a considerable 
disturbance in the gas disk of the target galaxy, we have 
found that it is rather ineffective in creating giant HI rings 
similar to that observed in Ho I. The HVC pushes the mat- 
ter mostly in the vertical direction so that the HI image 
of a nearly face-on galaxy does not show a considerable 
radial variation in the HI column density within the im- 
pact region. Indeed, Fig.ll7lshows the integrated HI image 
of the gas volume density distribution shown in Fig. EH 



An inclination angle of 5° is assumed. The formation of a 
ring-like structure is only seen in the late collision phase 
(t ~ 110 Myr), when the cloud has plunged through the 
gas disk and pushed a considerable fraction of the disk 
gas to an altitude z J> 1 kpc. Most of the gas consti- 
tuting this remnant retains the orbital momentum it had 
before the collision. Since the horizontal gravity is much 
less at z ^> 1 kpc than in the plane of the galaxy, the 
lifted gas starts spiraling away on ballistic orbits and it 
can be observed in nearly face-on galaxies as a dense HI 
ring surrounding the central depression (see the left lower 
panel in Fig I17|) . The maximum size of such an HI ring 
(diameter ~ 0.8 kpc) is at least twice smaller than that 
found in our numerical simulations of multiple SNe and 
hypernova explosions of the same total energy. Note that 
the diameter of the HI ring in Ho I is ~ 1.7 kpc. We find 
that most of the HI rings formed by the impact of HVCs 
have diameters less than 0.5 kpc. The HVCs of a smaller 
size and higher HI column density, but the same kinetic 
energy, tend to form smaller rings. The size of the HI ring 
depends weakly on the shape of the cloud. HVCs with the 
total kinetic energy below 2 x 10 53 do not form noticeable 
HI rings. HI rings with a diameter of > 0.8 kpc are never 
formed in our simulations for the considered upper limit 
of HVC kinetic energy, 2.5 x 10 54 ergs. 

A more dramatic difference between the HI rings 
formed by the vertical impact of HVCs and those created 
by multiple SNe is seen in the values of JC. As is obviously 
seen from Fig. HO 1000 consecutive SNe can create rings 
with K, as high as 100, whereas Fig. El indicates that the 
impact of a cloud of the same kinetic energy generates the 
ring with K. ~ 4.5. Note that /C in Holmberg I is ~ 15 — 20 
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Fig. 18. The pV-diagrams of the gas volume density dis- 
tribution shown in Fig. ^j] The pV-cuts are taken along 
the major axis of the projected galaxy. An inclination an- 
gle of 5° is assumed. 



(Ott et al. 2001). The same tendency is found for other 
values of released energy: multiple SNe produce rings with 
roughly 20-30 times higher values of K, than does the im- 
pact of a cloud of the same total energy. On the other 
hand, hypernova explosions form rings with K, roughly 
a factor of 2 higher than those created by the clouds of 
the same kinetic energy. The difference in /C between the 
rings formed by the three mechanisms considered above is 
smeared out as one considers higher inclination angles of 
i > 45°. 

The impact of HVCs should leave prominent features 
in the pV-diagrams, which could be used to distinguish 
between the rings formed by the impact of HVCs and 
those formed by multiple SN and hypernova explosions. 
Indeed, Fig. ^| shows the pV-cuts taken along the ma- 
jor axis of the galaxy experiencing the collision with a 
cloud of M H vc = 2.5 x 10 6 M© and w H vc = 200 km s" 1 
An inclination angle of 5° is adopted. The open bub- 
ble produced by the impinging cloud is clearly seen in 
the upper left panel of Fig. El The bubble starts filling 
out after t = 45 Myr. However, the cloud is still clearly 
visible in the pV-diagrams at that time as a dense ir- 
regular clump moving towards the observer. The veloc- 
ity of the cloud noticeably decreases during the passage 
from uhvc ~ 200 km s _1 at the onset of collision to 
vhvc ~ 20 km s _1 at the time when the cloud emerges on 
the other side of the gas disk. The pV-diagram at this time 
(t = 110 Myr) is dominated by two bright and compact 
clumps, which represent the cold dense gas compressed 
and pushed by the cloud to an altitude of z ~ 1 — 2 kpc. 



The depression in the HI flux seen in the pV-diagram near 
the galactic center is due to the HI ring formation dis- 
cussed above. Note that these two clumps are noticeably 
different in shape from those observed in the pV-diagrams 
of the shells formed by multiple SN and hypernova explo- 
sions (see Figs. I7I9|) . The second bubble formed by the 
remnant falling back onto the disk is clearly seen in the 
lower right panel of Fig. ^] 

It is obviously seen that neither the pV-diagrams nor 
the size and/or ring-to-center contrast /C obtained in the 
simulations of the vertical collision of HVCs with the 
galactic disk can account for what is actually observed 
in Ho I. We thus conclude that multiple SN explosions re- 
main the most plausible explanation for the origin of the 
giant HI ring morphology in Ho I. 

6. Discussion 

If SNe are the origin of HI rings surrounding HI holes, 
the presence of a remnant stellar cluster is expected. Such 
stellar clusters are not always observed within HI rings or 
holes. For instance, Rhode et al. (1999) found that several 
HI holes in the dwarf irregular galaxy Holmberg II lack 
any stellar cluster counterparts. However, in the case of 
Holmberg I, the B-band magnitude as well as the U — 
B color of the optical emission from the giant HI hole 
indicates that a stellar cluster of ~ 2 x 10 5 Mq may be 
present there (Vorobyov et al. (2004 ) . 

The high density clumps produced in our simulations 
due to the fragmentation of SNe-driven shells are usu- 
ally not self-gravitating. The masses of these clumps are 
slightly below the Jeans mass Mj = 5 c^d/ (4G), where 
d is the size of the clump. This expression is derived 
from the virial theorem assuming spherical clumps. This 
is consistent with a low Ha luminosity of Ho I: Lho = 
4.3 xlO 38 ergs s" 1 (Miller & Hodge 1996 ). This also implies 
that SNe alone have difficulty in producing self-gravitating 
clumps in dwarf irregular galaxies, and probably an exter- 
nal agent such as ram pressure or the collision of two or 
more SNe-driven shells is needed to induce star formation. 
In massive disk galaxies, radiative cooling is usually more 
efficient than in dwarf irregulars due to a higher metal- 
licity. Therefore, dense clumps are more likely to become 
Jeans unstable. On the other hand, the most energetic 
HVCs with kinetic energy > 1.0 x 10 54 ergs can indeed 
produce dense clumps that are Jeans unstable. However, 
such energies are rather exceptional and hence this mech- 
anism cannot account for most star formation in dwarf 
irregulars. 

There are a few assumptions inherent to the model 
that need further discussion. 

The interstellar medium distribution in Ho I. The solu- 
tion of the steady state momentum equation in Sect. |21 
produces a smooth initial gas distribution, which may be 
only a first order approximation to the more realistic in- 
homogeneous interstellar medium of dwarf irregulars. The 
influence of the deviations from a smooth ISM distribu- 
tion in the gas disk on the formation and evolution of 
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HI shells deserves a separate detailed study. We expect 
here that the effect of small scale inhomogeneities with 
sizes much smaller than that of the HI shell should be 
cancelled out as the shell forms. Conversely, large scale 
inhomogeneities with sizes comparable to that of the HI 
shell could influence its shape and evolution, making it 
easier/harder to break out of the disk. In particular, large 
scale inhomogeneities may reduce the smooth appearance 
of hypernova-driven shells discussed in Sect 14 .81 However, 
considering the size of the HI ring in Ho I (~ 1.7 kpc), 
we do not expect inhomogeneities of a such scale to be 
present in the initial gas distribution of this galaxy. 
The vertical scale height of the stellar disk and gas velocity 
dispersion. The actual value of the vertical scale height 
z s of the stellar disk is difficult to estimate, since Ho I 
has a low inclination angle. We have adopted a value of 
z s = 300 pc, which is typical for dwarf irregular galaxies. 
A moderate variation ±100 pc in the value of z s has a 
minor influence on the gas distribution and, hence, on our 
main results. The HI velocity dispersion of Ho I is mea- 
sured by Ott et al. I|2001l see their Fig. 10). A moderate 
value cthi ~ 9 km s _1 is found for most of Ho I disk, which 
is typical for gas-rich, quiet dwarf galaxies (^11 11999*1) and 
is used in our numerical simulations. A small portion of 
gas disk in the northwestern part of Ho I has, however, a 
slightly higher value of cthi ~ 12 km s _1 . Such a varia- 
tion in ohi is not expected to influence our main results. 
On the other hand, the measured cthi represents the cold 
component of the gas disk only. If a massive hot gas com- 
ponent is present in Ho I, the velocity dispersion cthi+hii 
of the mixture of hot and cold gas components may be 
significantly greater. An increased gas velocity dispersion 
makes the gas disk thicker. For instance, for an adopted 
value of cthi+hii = 20 km s _1 , the Gaussian scale height of 
the model gas disk h becomes roughly three times larger 
than that shown by the dotted line in Fig. However, 
h ^ 1 kpc is not expected in Ho I, because a SN-driven 
shell with radius r w 0.7 — 0.8 kpc would not blow out of 
the disk. Therefore, we would not observe such a high col- 
umn density contrast (JC « 15 — 20) between the HI ring 
and the central depression. Furthermore, Ho I is a quiet 
dwarf galaxy with an average star formation rate of only 
0.004 M Q yr _1 (Miller & Hodge UMty and we do not ex- 
pect a large amount of hot gas to be present there. Hence, 
we conclude that the integrated gas velocity dispersion of 
Ho I is not expected to deviate much from the adopted 
value of a = 9 km s -1 . 

Internal structure and geometry of HVCs. The shape 
of HVCs (i.e. spherical or cylindrical) is found to have 
a very minor influence on our results, which is agree- 
ment with many previous studies of HVC collisions with 
a gas disk (e.g. Comeron & Torra ITTJ92|) . The assump- 
tion of an initially homogeneous cloud is certainly an ide- 
alized case, since the observations show that HVCs are 
very clumpy and hierarchically structured (e.g. Wakker & 
Schwarz 1991). It may be that the remnant of the initially 
homogeneous cloud that emerges at the other side of the 
galactic plane at t rj 110 Myr in Fig.Elis a better example 



of an HVC prior to the collision. Inhomogeneous structure 
further reduces the ability of an HVC to produce large HI 
rings. 

7. Conclusions 

We have numerically studied the formation of a giant 
HI ring as is observed in the Holmberg I dwarf irregu- 
lar galaxy (Ott et al. 12001(1 . The following three energetic 
mechanisms with the total released energy of <J 10 54 ergs 
have been considered: multiple SNe, a hypernova explosion 
associated with a gamma ray burst, and the vertical im- 
pact of a HVC. The shells created by those mechanisms 
are seen in the integrated HI image of a nearly face-on 
galaxy as the HI rings of various sizes and ring-to-center 
contrasts JC in the HI column density. We find the fol- 
lowing differences among the HI rings formed by those 
mechanisms: 

• Multiple SNe form rings with the highest ring-to-center 
contrast JC <; 100. A single hypernova and an HVC form 
the rings with a much smaller values of JC <J 10 and K, <^ 5, 
respectively. This is due to the difference in the dynamics 
of the shells produced by multiple SNe and those created 
by a hypernova: multiple SNe are much more effective in 
producing blowout due to a combined action of Kelvin- 
Helmholtz and Rayleigh- Taylor instabilities than are sin- 
gle hypernovae of the same total energy. The vertical im- 
pact of an HVC is found to be ineffective in evacuating 
the gas in the radial direction and hence in creating HI 
rings with a high ring-to-center contrast. It is important 
to note that the difference in K. found among three en- 
ergetic mechanisms smears out at an inclination angle of 
i > 45° due to the projection effect. 

• Both multiple SNe and a hypernova can account for the 
observed size of the HI ring in Ho I (diameter ~ 1.7 kpc). 
The maximum size of the ring that an HVC is found to 
form in our numerical simulations is ~ 0.8 kpc. We thus 
conclude that only multiple SNe can reproduce both the 
size (diameter ~ 1.7 kpc) and the ring-to- center contrast 
in the HI column density (JC ~ 15 — 20) of the HI ring in 
Ho I. 

We construct the model position-velocity diagrams 
and find that they can be used to distinguish between 
HI rings produced by different energetic mechanisms. The 
model pV-diagrams of the blownout shell produced by 
multiple SNe are similar to those observed in Ho I and 
show two symmetrically displaced, elongated blobs rep- 
resenting the walls of the shell. On the other hand, the 
model pV-diagrams of the shell created by a hypernova 
explosion show a characteristic elliptical ring-like struc- 
ture, which is not seen in the observed pV-diagrams of 
Ho I. The impact of HVCs leaves prominent features in 
the model pV-diagrams that are also not found in the ob- 
served pV-diagrams of Ho I. Those features include dense 
irregular clumps and a void representing the impinging 
cloud and the bubble formed by the impact. We thus con- 
clude that multiple SNe appear to be the most plausible 
energetic mechanism of the HI ring formation in Ho I. 
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Our numerical simulations indicate that the appear- 
ance of the SN-drivcn shell in dwarf irregulars may depend 
on the inclination of the galaxy. In nearly face-on galaxies 
the blownout shell appears in the integrated HI image as 
an HI ring with a deep central minimum. There is little 
azimuthal variation in the HI column density around the 
ring. The integrated HI image of the blownout shell in 
considerably inclined galaxies [i > 45°) is instead dom- 
inated by two kidney-shaped enhancements and shows a 
mild central minimum, which is similar to what is found 
in Sectans A and M81 dwA (Skillman et al.[I988| Stewart 
1998). In nearly edge-on galaxies the integrated HI image 
of the blownout shell resembles a 'dumbbell'. 

The offplane explosion of a hypernova creates the egg- 
like shell that can preserve its shape for a considerably 
longer time than the shell created by multiple SNe of 
the same total energy. This may explain the smooth ap- 
pearance of the highly elongated Ha shell observed by 
Dennison et al. l|1997f) in the W4 HII region of the Perseus 
arm. 
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